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Abstract 

CD 

$—1 ' Inelastic neutron scattering is applied to study the role of magnetism in stabilizing the charge ordered 



State in /?i/3Sr2/3Fe03 (RSFO) (R = La, Pr, and Nd). The ratio of the ferromagnetic exchange energy 
(Jp) and antifen^omagnetic exchange energy (Jaf), \Jf/Jaf\, is a key indicator of the stability of the 



'^ ■ charge ordered and antiferromagnetic ordered state. This ratio is obtained from the spin wave spectrum 

O ■ by inelastic neutron scattering and is sufficiently large to suggest that the magnetic exchange energy 

alone can stabilize the charge ordered state in Lai/3Sr2/3Fe03 and Pri/3Sr2/3Fe03. The exchange ra- 
'i> I tio decreases from Lai/3Sr2/3Fe03 to Ndi/3Sr2/3Fe03 which indicates a gradual destabilization of the 

O ■ magnetic exchange mechanism for charge ordering in correspondence with the observed reduction in the 

^i . ordering temperature. 

o. 

PACS numbers: 



X 



I) Introduction 

The charge-ordering (CO) transition is often encountered in complex transition-metal oxides 
(TMO) and has been the focus of intense inquiry and debate in condensed matter science in the 
past years. The metal-insulator transition that occurs as the temperature decreases across the 
CO transition temperature, Tco, is associated with a change from an itinerant electronic state to 
a more localized state. The CO state plays an important role in various systems, including the 
superconducting cuprates and the proximity of the superconducting state to a spin/charge stripe 
ordered state,[l] colossal-magnetoresistive manganites, where the CO states compete with fer- 
romagnetic metallic states, [2] and layered nickelates, which form a small polaron lattice. [3] 
Therefore, understanding the causes and implications of CO phenomena is significantly impor- 
tant. The CO state is also often closely associated with magnetic and orbital ordering, and it is 
widely recognized that CO can arise from a variety of competing interactions, most importantly 
the intersite Coulomb interaction, magnetic exchange energy, and electron-phonon interactions, 
all of which are strongly dependent on the valence states of neighboring metal ions. 

The Sr-doped rare earth ferrite i?i/3Sr2/3Fe03 (7?SF0) is an interesting example of a CO 
system where magnetic exchange energy is thought to play a crucial, if not dominant, role in 
the stability of the CO state. i?SFO is a perovskite based crystal where the Fe ion adopts a formal 
fractional valence of 3.61+. Below Tqo^ it has been proposed that charge disproportionation 
occurs according to 3Fe'^^^"*" ^ 2Fe'^''' -i- Fe^+, with the different iron valences ordering in 
planes containing a repeating arrangement of 3-I-, 3-I-, 5-i- ions perpendicular to the body diagonal 
[lll]c.[4, 5] The CO occurs simultaneously with antiferromagnetic (AF) order. Recently, we 
used inelastic neutron scattering (INS) measurements of the spin wave spectrum to demonstrate 
the plausibility that the magnetic exchange energy is the dominant interaction giving rise to 
CO in Lai/3Sr2/3Fe03 (LSF0).[6] Our results show that the observed CO ground state can be 
stabilized by large ferromagnetic (F) exchange {J p) occurring between the nearest-neighbor 
(AW) Fe^"*"-Fe^"^ pairs. However, this conclusion was made based on the assumption that the 
intersite Coulomb interaction (i.e. the Madelung energy) is suppressed by strong electronic 
screening enabled by the small charge transfer (CT) gap (several to lO's of meV) observed in 
the system. [6] Since we cannot directly determine the contribution of the Coulomb interaction, 
one possible way to verify the dominance of the magnetic exchange mechanism is to measure 
the magnetic exchange energies in other i?SFO (for example, 7? = Pr and Nd) compounds. For 



the smaller Pr^+ and Nd^+ ions, increased lattice distortions lead to a larger charge-transfer gap 
and a narrower electronic bandwidth. This should lead to reduced screening and increase the 
stability of the CO state due to a greater influence of the Coulomb interaction. However, r^v 
and Tco are known to be suppressed by smaller i?'^+. [5] In 1998, T. Mizokawa and A. Fujimori 
proposed that the \JfIJaf\ exchange ratio is a good indicator of the propagation direction of 
CO ordering in the limit where magnetic energy is dominant. If the ratio is larger than 1, the 
charges will be ordered along [111]; if the ratio is less than 1, the charges will be ordered along 
[100]. [7] Thus, if neutron scattering measurements can indicate a weakening of the magnetic 
exchange ratio, this would give additional support to the magnetic mechanism for CO. 

To verify this hypothesis, we study the spin wave spectrum in 7?SF0 with different R^^ 
ions. Based on INS measurement of the powder sample in principal, the magnetic exchange 
energies, Jp and J af, can be obtained,[8] and the related exchange ratio could be calculated. 
The corresponding contribution of the magnetic energy to the CO state could then be considered. 
However, the magnetic spectra of PSFO and NSFO are not as simple as that of previously 
studied LSFO system,[6] where only Fe moments contribute to the neutron magnetic cross- 
section, as the Pr^+ and Nd^+ ions themselves possess magnetic moments. The neutron intensity 
from crystal electronic field (CEF) excitations of the magnetic Pr^+ and Nd^+ ions makes the 
extraction of the Fe spin wave spectrum more complicated. To better account for the rare earth 
CEF excitations, we also investigate the parent compounds, RVO (R = Pr and Nd). All of 
the i?FO parent compounds are insulators with G-type antiferromagnetic ordering of Fe atoms 
occurring at high temperatures (T^ ^ 700 K). The rare earth ions in RFO are expected to have 
similar CEF spectra to RSFO once the simple dilution of the rare earth site by Sr is taken into 
account. After accounting for the CEF intensities of the R^^ ions in the total magnetic cross- 
section of RSFO, the INS data may be compared to calculations of the spin wave spectra and 
their corresponding cross-sections using a Heisenberg model. We show that the ratio | JfUaf I 
is found to decrease with smaller i?^+ which can account for the reduction in Tj^ by the magnetic 
mechanism. 

II) Experiment 

Polycrystalline i^FeOs {RFO) and RSFO (R = La, Pr, and Nd) were prepared by a conven- 
tional solid-state reaction method. Stoichiometric amounts of La203, or Nd203, or PreOn, 



SrCOa, and Fe203 were mixed by grinding with mortar and pestle. The mixtures were trans- 
ferred to an AI2O3 crucible and calcined several times in air at temperatures of 1100 °C and 
1200 °C respectively for 24 hours. Then, the press-formed pellets were sintered in air at 1250°C 
and 1350°C for 30 hours, respectively. As the ionic size decreases from La to Nd, the RSFO 
compounds tend to be more oxygen deficient. [9] PSFO and NSFO were annealed under oxygen 
pressure (10 bar) at 600°C for 72 hours. Room temperature powder X-ray diffraction (XRD) 
patterns were performed on a Rigaku Miniflex X-ray diffractometer with Cu Ka radiation to 
confirm phase purity. No impurities were observed. The oxidation state of iron was determined 
by iodometric titration and is listed in Table I. 
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FIG. 1: (color online) Magnetization of 7?SF0 {R = La(a), Pr(b) and Nd(c)) as determined by ZFC 
SQUID measurements; (d) powder neutron diffraction of PSFO from 300 K to 15 K as a function of 
d-spacing. The arrows show the positions of the magnetic Bragg (red) and charged-ordered superlattice 
(black) peaks. 



The CO and AF transition temperatures for the Sr-doped samples were determined by neu- 
tron powder diffraction (NPD), using the High-Intensity Powder Diffractometer (HIPD) at the 
Lujan Center at Los Alamos National Laboratory, and zero field cooled (ZFC) magnetization 
measurements using a superconducting quantum interference device (SQUID) magnetometer. 
These characterization data are shown in Fig. 1 . The RSFO systems were found to have iden- 
tical magnetic propagation vectors and charge-order propagation vectors (inferred via the de- 
velopment of structural superlattice peaks) below T^- Fig. 1(d) shows temperature dependent 
neutron diffraction data for PSFO and indicates the charge-order propagation vectors at (n/3, 
n/3, n/3), and magnetic propagation vectors at (n/6, n/6, n/6).[10] The lattice constants deter- 
mined from refinement of 300 K NPD and magnetic transition temperatures are listed in Table 

I. The bond lengths and bond angles determined from NPD patterns at 300 K are listed in Table 

II. And the geometric tolerance factor (t) is expressed as, 

<R-0> ^,, 



V2<Fe-0>' 
where <R-0> and <Fe-0> are the average bond lengthes of R-0 and Fe-0. 

INS measurements were performed on the Pharos spectrometer at the Lujan Center of Los 

Alamos National Laboratory and the ARCS spectrometer at the Spallation Neutron Source 

(SNS) of Oak Ridge National Laboratory. Both instruments are direct geometry time-of- flight 

(TOF) spectrometers and measure the scatterring intensity over a wide range of energy transfers 

(hcu) and scattering angles between 1°-140°, thereby allowing determination of a large swath of 

the scattering intensity, 5(Q, u), as a function of momentum transfer (hQ) and energy transfer 

(hio). 

On ARCS, powders (~14 g) of R¥0 and i?SFO (R = La, Pr, and Nd) were packed in 5 
aluminum foil sachets. The sachets were placed in an aluminum can filled with He exchange 
gas whose size was approximately 4.5 cm x 6.5 cm x 0.5 cm. INS spectra were measured 
with an incident energy (Ei) of 180 meV. On Pharos, ~ 50 g of 7?F0 and RSFO {R = La, and 
Nd) were packed in a flat aluminum can (6 cm x 6 cm x 0.5 cm), and EiS were 120 meV and 
160 meV. The face of the sample can was oriented at 135° to the incident neutron beam for 
both instruments in a transmission geometry. To achieve adequate statistics, the sample was 
measured for approximately ~ 24 hours on Pharos, and ~ 5 hours on ARCS. Empty sample 
can measurements were also performed and subtracted from the data presented. 



TABLE I: Lattice and magnetic parameters of RFO and RSFO (R = La, Pr and Nd) as determined by 
X-ray/neutron scattering measurements at 300 K, SQUID measurements(riv), and iodometric titration 
(oxygen deficiency). 





LaFeOs 


PrFeOa 


NdFeOg 


Space group 


Prima 


Pnma 


Pnma 


Lattice constant(A) 
a 
b 
c 


5.56 ± 0.01 
7.85 ± 0.02 
5.56 ± 0.01 


5.57 ± 0.01 
7.79 ± 0.02 
5.48 ±0.01 


5.59 ± 0.01 
7.76 ± 0.02 
5.45 ± 0.01 


r7v(K)[ii] 


738 


707 


693 




Lai/gSrs/aFeOa 


Pri/sSra/sFeOg 


Ndi/gSra/gFeOg 


Space group 


/?3c 


/?3c 


/?3c 


Lattice constant(A) 

a-b 
c 


5.48 ± 0.01 
13.41 ± 0.03 


5.48 ± 0.01 
13.37 ± 0.03 


5.47 ± 0.01 
13.34 ± 0.03 


riv(K) 


210 ±2.0 


191 ±2.0 


189 ± 2.0 


Oxygen stoichiometry 


2.94 ± 0.03 


2.97 ± 0.03 


2.97 ± 0.03 



III) Results and Discussion 



The unpolarized inelastic neutron scattering cross-section contains contributions from both 
magnetic and phonon scattering. In order to isolate the spin wave spectrum, the magnetic 
scattering must first be separated from the phonon scattering. This is accomplished by using 
the fact that the magnetic scattering intensity decreases with Q (or 26) due to the magnetic form 
factor, while phonon scattering intensity increases proportional to Q^. 

The INS data of NFO taken on ARCS at T = 10 K with £,; = 180 meV are shown in Fig. 
2(a) and used as an example of the data treatment. A similar analysis was performed for LFO 
on Pharos as outlined in Ref. [6]. The data summed over the high angle range of 26 = 75- 
95° contain primarily phonon scattering. Fig. 2(b), while the data within the low angle range 
of 10-30° contain scattering from both phonons and spin waves arising from the G-type AFM 
order of Fe^+ and the CEF of Nd'^+, Fig. 2(c). The magnetic scattering in NFO was isolated 



TABLE II: The structural properties of RFO and RSFO as determined by NPD at 300 K. 

LaFeOs PrFeOs NdFeOs 
Bond length(A) 

<R-0> 2.694 ± 0.008 2.627 ± 0.008 2.609 ± 0.008 

Fe-O(l) 2.002 ±0.004 2.004 ± 0.004 2.005 ± 0.004 

Fe - 0(2) 2.004 ± 0.004 2.006 ± 0.004 2.007 ± 0.004 

Fe - 0(2) 2.005 ± 0.004 2.015 ± 0.004 2.017 ± 0.004 

Bond angle 

ZFe-0(l)-Fe 157.6° ± 0.3° 153.3° ± 0.3° 151.2° ±0.3° 

ZFe-0(2)-Fe 157.5° ± 0.3° 152.4° ± 0.3° 151.4° ±0.3° 

Geometric tolerance factor 

t 0.951 ± 0.001 0.925 ± 0.001 0.918 ± 0.001 

Lai/gSra/sFeOa Pri/gSrs/sFeOs Ndi/gSra/sFeOa 
Bond length(A) 

R - O 2.7413 ± 0.006 2.7377 ± 0.006 2.7329 ± 0.006 

Fe-0 1.940 ±0.004 1.941 ±0.004 1.939 ±0.004 
Bond angle 

Z Fe - O - Fe 173.2° ± 0.4° 170.5° ± 0.3° 169.3° ± 0.3° 
Geometric tolerance factor 

t 0.999 ±0.002 0.997 ±0.002 0.996 ± 0.002 

by subtracting the high angle phonon data from low angle data after scaling the high angle data 
by a constant factor. A comparison of the scaled high angle data to the low angle data is shown 
in Fig. 2(c) and the phonon subtracted data is shown in Fig. 2(d). In order to compare the two 
instruments used for the INS measurements, the magnetic spectrum of NFO measured on the 
Pharos is overplotted in Fig. 2(d). The spectrum agrees with each other very well. 
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FIG. 2: (color online) (a) Inelastic neutron scattering intensity (Ei =180 meV) of NFO (color scale) 
versus scattering angle and energy transfer at T = 10 K as measured at ARCS. Horizontal white lines 
delineate regions where phonon and magnetic scattering are isolated; (b) neutron intensity summed over 
the angular range from 75-95° originating from phonons (red dots); (c) neutron intensity summed over 
the low angular range from 10-30° (blue dots) and phonon background scaled from the high angle sum; 
(d) the isolated magnetic scattering of ARCS data (red dots) and Pharos data(black line). 

A) Magnetic spectra of /?FO 
a) Fe Spin Waves 



We now discuss the analysis of collective Fe spin waves below T^ in the parent RFO com- 
pounds. In RFO, NN Fe'^+(3J^) spins are coupled by strong AF superexchange interactions 



{Jaf < 0). According to the single-crystal INS studies of TmFe03,[12] the spin waves can be 
approximated using a Heisenberg model Hamiltonian with only isotropic AW exchange interac- 
tion. 



H 



~Jaf / ^ Si ■ Sj , 



(2) 



<i,j> 



where Sj and Sj represent the spin vectors of the ith andjth iron atoms that are NNs. 

The INS information obtained from the polycrystalline samples is related to the spin wave 
density-of-states (SWDOS) via a powder-averaging over all crystallographic directions. [8] For 
the G-type LFO spin waves, the SWDOS consists of a single sharp peak at an energy of6\JAF \ 
S^~^ and S'^"'" is the spin magnitude of Fe'^"'^ ion. Assuming S^^ = 5/2, a value of Jaf = -4.9 meV 
can be determined from the position of this single peak at 73 meV, shown in Fig. 3. The value 
of Jaf is also listed in Table III. 
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FIG. 3: (color online) The experimental angle-averaged magnetic INS data for RFO (R = La, Pr and 
Nd) at lOK (dots). Heisenberg model calculation of the powder-averaged Heisenberg spin waves cross- 
section (Heisen) for G-type magnetic order are shown as solid cyan Unes. Additional peaks in the spectra 
for R = Pr and Nd and associated solid pink lines are the fits of CEF excitations. 



Figure 3 shows the phonon subtracted magnetic data for LFO, PFO, and NFO found via the 
procedure described in Fig. 2. As La"'"'^ has no f-electrons, there are no CEF excitations existing 
in LFO. While PFO and NFO have similar G-type magnetic structure as LFO, it is clear that 



the neutron spectra of PFO and NFO contain addiontal magnetic excitations due to R^^ CEP 
excitations which will be discussed later. In order to determine which of the observed peaks 
are due to Pe spin waves and which are CEPs, we can use mean-field theory and knowledge 
of Tn to estimate Jaf exchange coupling between Pe ions for R = Nd and Pr. In mean-field 
theory, [13] 



SknTM = N\Jaf\ VSiS + 1) , (3) 

where N = 6 and S is the spin angular momentum, 5/2. [8] With this assumption, the magnetic 
exchange energies are expected to weaken slightly upon going from LPO to NPO due to the 
decrease of T^. However, it is well-known that T^- is overestimated in mean field theory and 
therefore the mean-field estimate is smaller than that observed by neutron scattering. Table 
III lists the AP exchange energy based on Eq. (3). A more accurate value of the exchange 
can then be obtained by fitting the SWDOS peak in the INS data by Heisenberg model, Eq. 
(2). The results are listed Table III and fits to the SWDOS are shown in Pig. 3. According to 
superexchange theory, crystalline distortions caused by the smaller rare earth ions result in a 
bending of the Pe-O-Pe bond angle that weakens Jaf, Table II. [14] 

TABLE III: Magnetic exchange energy of 7?F0 (R = La, Pr and Nd). 



Jaf (meV) 


LaFeOs 


PrFeOs 


NdFeOa 


Mean Field 


-3.74 


-3.58 


-3.51 


Heisenberg model 
(based on LFO) 


-4.9 


-4.69 


-4.60 


Heisenberg model (fit) 


-4.9 


-4.55 


-4.45 



b) Crystal Electric Fields of Pr3+ and Nd^+ 

Mean field theory helps us locate the characteristic energy of spin-wave excitations from the 
Pe"^"*" ions. The remaining excitations in the phonon subtracted data should be from CEP of i?^"*" 
ions. The cross section for CEP excitations can be written as, [15] 

^'"^ oc[gjF{Q)re-'^^ScFF{Q,co), (4) 



d^du k 



10 



where Q = k' - kis the scattering vector, and fiio is the energy transfer, gj is Lande factor of the 
R^^ ion, F(Q) is the R^^ magnetic form factor, k' and k are values of the initial and final neutron 
wavevectors, and e~^'^ is the R^~^ Debye-Waller factor. ScEpiQ, ^) is the response function 
of the system which is determined by the temperature, CEF level splitting, CEF eigenstates and 
corresponding CEF matrix elements. 

Scef{Q,oj) = J2p,\{t\J^\ j) |2 SiE, - E^ - hu) , (5) 

where | / ) and | j ) are the initial and final CEF eigenstates of the system with level energies 
Ei and Ej. J± is the component of the total angular momentum operator perpendicular to the 
scattering vector; pi is the thermal population factor of the initial state. Observable excitations 
occur between levels which have non-zero matrix elements. 

For rare earth ions, the spin-orbit coupling is usually stronger than the CEF potential, and the 
total angular momentum J = L -i- S, remains a good quantum number. Therefore, the magnetic 
form factor in Eq. (4) is given by, [15, 16] 

where (jo) and (72) are Q-dependent functions whose values are tabulated. [17] 

b.l) NdFeOa 

The ground state of the free Nd^"*" ion has 3 unpaired /-electrons, and the Russel-Saunders 
term symbol is '^Iq/2 with a degeneracy of 2 7 -1- 1 = 10. At a point of orthorhombic symme- 
try in the distorted perovskite cell, the ground state multiplet splits into (2 J + l)/2 = 5 CEF 
doublets. [16] At 10 K, we observe four CEF transitions at ~ 9, 21, 46 and 60 meV as shown 
in Fig. 3. These are consistent with previous work,[18, 19] and are associated with excitations 
from the CEF ground state to each of the four excited states. Due to the possible overlap of 
CEF excitations with spin waves, the phonon scattering contribution and multiple scattering, 
the separation of the CEF contribution was done by examining both the u- and Q- dependence 
of the total cross-section. The positions and integrated intensities of the Nd'^+ CEFs in NFO at 
10 K were determined by. 



11 



S{Q, u) — SmagiQ, ^) + SphononiQ, ^) + ScEf{Q, ^) + Smulti{Q, ^) , (7) 

where SmagiQ, ^) is the polycrystalline averaged spin wave scattering of Fe^+ ions, SphononiQ, 
uj) is the polycrystalline-averaged phonon background, and Scef (Q, ^) is CEF scattering from 
Nd^"*" ions. 

For simplicity, we treat the phonon scattering from a powder sample in the incoherent ap- 
proximation. In this approximation, the one-phonon scattering is proportional to the phonon 
density-of-states (DOS) and can be expressed as, 

Sinc,±l-phonon = —Z{flu){n + 1) , (8) 

where X{h u) is the sum of weighted partial phonon DOS, Zj(^ uj), 

Z{nw) = V ^e-^^^Zi{hu) , (9) 



^2M, 



and (n+1) is the Bose population factor, 

{n + l) = ^[coth{^nwP) + l], (10) 

Constant energy scan obtained from averaging over different energy transfer ranges are 
shown in Fig. 4. Figure 4(a) is averaged over an energy range from 10.5-20.5 meV and 
described here as an example. The multiphonon contribution and other background contri- 
butions (multi.) are treated as a constant background (horizontal line). The single phonon 
(ph.) cross-section is determined by assuming an incoherent quadratic Q-dependence (parabola 
curve) given by Eq. (8). We note that the peaks seen at high-Q (above 5 A^^) arise from coher- 
ent phonon scattering that is not included in the analysis. The spin wave scattering intensities 
(Heisen) of Fe^"^ are calculated from a Heisenberg model (zigzag curve). Note that the sharp 
peaks in the spin wave contribution arise from the coherent scattering, which is included in the 
Heisenberg model calculations of the powder averaged cross-section. Especially at low ener- 
gies, the sharp inelastic peaks are coincident with the position of magnetic Bragg peaks. The 
remaining signal is associated with CEF scattering, which follows the magnetic form factor. 
The energy ranges where there is a large difference between total fitting line (multiphonon -i- 
phonon -i- spin wave) and the data signals the presence of a CEF excitation. The dotted line in 
Fig. 4 is the fitted CEF excitations. 
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FIG. 4: (color online) The Q-dependence of the neutron scattering data averaged over different energy 
transfer ranges in NFO at T = 10 K. (a) 10.5-20.5 meV, (b) 30.5-40.5 meV, (c) 40.5-50.5 meV, and (d) 
50.5-60.5 meV. The black dots are the experimental data. The blue line is an estimate of the incoherent 
phonon background {ph.) plus multiple scattering (multi.). The green line is the calculation of the poly- 
crystalline averaged spin wave scattering (Heisen) plus background using the parameters in the text. The 
red line is sum of the calculation of CEF, the polycrystalline averaged spin wave scattering, background, 
and multiple scattering. The vertical brown dash line in (b) - (d) is the kinematic limit for neutron data. 
The dotted line is the fitted CEF excitations. 

The measured and calculated energies and transition intensities of the CEF levels are listed in 
Table IV and plotted in Figs. 3 and 4. The CEF excitations of NFO can be found via the single- 
particle crystal field theory, and is discussed in Refs. [18, 19]. The calculation and measurement 
of the CEF transitions agree with each other very well. If we normalize the intensities of the 
CEF transitions to the intensity of the excitation from the ground state to the first excited state, 
which are 9.4 and 10.1 meV for the measurement and calculation respectively, the measured 
intensities of the second ( ~ 21.2 meV) and third excited states (~ 45.7 meV) agree with the 
calculations based on estimates of the corresponding matrix elements. A comparison of the 
intensity of the 59.9 meV is more difficult due to the proximity to the magnetic signal from 
Fe'^+ spin waves. 
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TABLE IV: The measured and calculated CEF transition energies of Nd^+ in NFO at 10 K.[18, 19] 



Energy levels 


fyjj{ 


meV) 


Integrated Intensities 


calculated 


measured 


calculated 


measured 





0.0 


0.0 






1 


10.1 


9.4 ±0.1 


100.0 


100.0 ±3.1 


2 


22.4 


21.2 ±0.1 


29.3 


30.9 ±1.1 


3 


44.7 


45.7 ±0.1 


46.1 


41.8 ±1.6 


4 


60.8 


59.9 ±0.1 


7.0 


17.8 ± 1.3 



b.2) PrFeOs 

We next examine the CEF transitions in PFO. Although many experimental and theoretical 
studies on Pr'^+ in perovskite oxides have been performed, [20, 21] we know of no complete set 
of experimental data for CEF levels of Pr'^"'" ions in the perovskite structure ABO3. There are 
2 unpaired /-electrons for Pr'^+ and the electronic term is ^H^. Since the structure of PFO is 
orthorhombically distorted and / = 4, the ground state multiplet splits into 27 + 1 = 9 singlets. 
The dipole allowed transitions from the ground state are listed in Table V. Similar to the analysis 
of NFO, the magnetic form factor of Pr^+ in PFO is taken from the literature[17] and used 
with the previously observed transitions from Ref. [20] in Table V to establish the fraction of 
different contributions to the total cross-section. Fig. 5 shows the result of the process. The 
data agree well with the model calculation over a range of wave vector and energy transfers. 

The CEF levels in PFO were compared to those of Pr'^^ in PrGaOs (PGO) because it had 
a similar Pr'^+ CEF scattering as PFO and was the only theoretical calculation on CEF of Pr^+ 
in perovskite oxides with the space group Pnma.[20, 21] There are discrepancies between our 
measurements of CEF excitations in PFO and the calculations for PGO. First, the CEF transition 
predicted to appear around ~ 67 meV was not observed in our measurement. This is likely 
due to the proximity of this transition to the Fe^+ spin wave band. The proposed observation 
of this CEF transition was mentioned in regards to INS measurements of PFO and PGO. [20, 
21]. However, this analysis did not report the spin-wave excitation of Fe'^+ ions, [21] and the 
magnetic peak claimed to be a CEF around ~ 60 meV is much more likely to be Fe'^+ spin 
waves scattering. An examination of Fig. 5(e) shows that the 67 meV CEF transition may 
be either too weak to observe (as compared to the Fe^+ spin waves) or shifted in energy as 
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FIG. 5: (color online) The Q-dependence of the neutron scattering data for different energy transfer 
ranges (Q-cuts) in PFO at T = 10 K. The brown dash line in some different energy transfer ranges ((c), 
(e), and (f)) indicates the usable low-angle data limit of the experiment. The dotted line is the fitted CEF 
excitations. 

compared to PGO. We also observe a weak excitation peak around ~ 100 meV, Fig. 3. This is 
presumably a CEF excitation, however the matrix element is predicted to be zero. [20] 

B) Magnetic excitations in RSFO 

We now take our knowledge of the different scattering cross-sections in the parent RFO 
compounds and use it in an attempt to isolate the scattering arising from the Fe"'^ spin waves 
in the doped RSFO compounds. Figure 6(a) shows the extracted low-angle magnetic intensity 
of 7?SF0 (R = La, Nd, and Pr). In general, we find magnetic signals up to a maximum energy 
of 120 meV. The high energy portion of the magnetic excitation spectrum between 90 ~ 120 
meV is very similar for each compound. In LSFO, this high energy band is associated with spin 
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TABLE V: The measured and calculated CEF transition energies of Pr^+ in PFO at 10 K. 



Energy levels 



FrFeOs 

(measured) 



PrGaOs Integrated Intensities 

(calculated) (present work) 




1 

2 
3 
4 
5 
6 
7 
8 



Ref. [20] present work Ref. [21] 



0.0 


0.0 


0.0 


2.0 ±0.1 


-- 


5.6 


14.7 ± 0.4 


15.2 ±0.1 


16.0 


23.2 ± 0.5 


24.7 ±0.1 


23.4 


36.0 ± 1.0 


36.6 ± 0.5 


32.9 


58.0 ± 2.0 


-- 


67.4 


-- 


-- 


69.3 


-- 


80.1 ±0.5 


89.5 


__ 


97.9 ± 0.3 


113.1 



0.0 

100.0 ±3.1 

41.1 ± 1.6 

82.2 ± 2.3 



4.2 ± 0.8 
4.9 ± 0.5 



waves propagating in ferromagnetic channels in the magnetic lattice. The low energy features 
are more difficult to compare due to the presence of R^^ CEF excitations. 

a) Crystal Electric Fields of Pr3+ and Nd3+ 

The CEF information obtained from the parent RFO compounds (R = Pr and Nd) can be 
used as a guide to estimate the CEF contribution in Sr-doped RSFO. Because the ionic radius of 
Sr^+ is close to the radii of the R^^ ion, we assume that the average structural environment of 
the i?^+-site does not change significantly in RSFO with Sr^+ doping. Thus, the CEF of i?^+ in 
RSFO would be similar to that of 7?F0. Since 2/3 of R^^ has been substituted by nonmagnetic 
Sr^+ ions, the integrated intensity of the CEF excitations in RSFO should be 1/3 of that found 
for RFO. However, the effect of disorder and other lattice distortions arising from the Sr^+ 
substitution may shift the position and broaden the width of CEF excitations. 

The CEF scattering intensities of RSFO were fitted with the same method as used for the 
RFO samples. The intensities were initially constrained to be one-third of the corresponding 
transition in RFO and then allowed to vary in the final fits. The fitted results of R^^ CEFs in 
the compounds RSFO are compared with the data of i?FO in Fig. 6(b) and (c). The integrated 
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FIG. 6: (color online) (a) Magnetic inelastic neutron scattering intensity of RSFO {R = La, Pr, and Nd) 
versus energy transfer at T = 10 K and Ei= 180 meV. The calculated CEF intensities of (b) PFO (black 
line) and PSFO (red line); (c) NFO (black line) and NSFO (red line), (d) The spin wave scattering from 
Fe"^"^ ions in RSFO (R = La (black), Pr (red), and Nd (blue))found via a difference of the total magnetic 
signal (a) minus the fitted CEF contribution of (b) and (c) respectively. 

intensities of R^^ CFE states in these compounds are listed in Table VI (Pr^+) and Table Vn 
(Nd"^"^), respectively. The transition energy between different energy levels is defined as Eji, 
which ith and j'th are the numbers of the energy levels. 

The integrated intensity ratio of the CEF excitations of Nd'^"'" in NFO is consistently 3 times 
that of the CEFs measured in NSFO as shown in Table VIL For PSFO, the ratio of Pr'^^ CEF 
intensities deviates from 3 for the £'20 transition, which overlaps the strong elastic peak. 
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TABLE VI: The transition energies and integrated intensities (area) of Pr'^+ CEF transitions in PFO and 
PSFO at 10 K. The last column is a ratio of the integrated intensities of the related excitations. 



Ejo 


PFO 


PSFO 


PFO/PSFO 


fvjj (meV) area 


hw (meV) area 


Eio 


__ 


__ 


-- 


E20 


15.2 ±0.1 1.29 ±0.04 


13.9 ±0.3 0.68 ±0.01 


1.9 ±0.1 


E30 


24.6 ±0.1 0.53 ±0.02 


23.3 ±0.2 0.19 ±0.01 


2.9 ± 0.2 


E40 


38.6 ±0.5 1.06 ±0.03 


34.1 ±0.1 0.36 ±0.01 


2.9 ±0.1 


E50 


__ 


__ 


-- 


Eeo 


__ 


__ 


-- 


E70 


88.1 ±0.5 0.054 ±0.01 


85.9 ± 0.4 0.018 ± 0.003 


3.0 ± 0.4 


Ego 


97.9 ± 0.3 0.063 ± 0.07 96.9 ± 0.4 0.021 ± 0.002 


3.0 ± 0.4 



TABLE VII: The transition energies and integrated intensities (area) of Nd^"*" CEF transitions in NFO 
and NSFO at 10 K. The last column is a ratio of the integrate intensities of these excitations. 



EjO 


NFO 


NSFO 


NFO/NSFO 


E (meV) area 


E (meV) area 


Eio 


9.1 ±0.1 1.8 ±0.1 


9.1 ±0.4 0.6 ±0.1 


3.0 ± 0.3 


E20 


20.8 ±0.1 0.6 ±0.1 


19.2 ±0.2 0.2 ±0.1 


3.0 ± 0.4 


E30 


45.2 ±0.1 0.9 ±0.1 


46.9 ±0.4 0.3 ±0.1 


3.0 ± 0.5 


E40 


59.5 ±0.2 0.5 ±0.1 


60.9 ±0.4 0.2 ±0.1 


3.1 ±0.8 



b) Fe Spin Wave Excitations 

After subtracting the CEF intensities for the 7?SF0 compounds, the spin wave excitations of 
Fe ions in all three 7?SF0 compounds were isolated and are compared in Fig. 6(d). The spin 
wave spectra in all three compounds agree with each other for energies greater than 85 meV. 
At all other energies, the spin wave spectra of LSFO and PSFO remain similar: there are two 
energy bands with a gap between 60-80 meV. In NSFO, the low energy spectral weight appears 
to move to higher energies and fills this portion of the spectrum. The intense peak around 20-30 
meV in NSFO is a signal which also exists in both LSFO and PSFO and is likely an artifact due 
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to inaccuracies in subtracting the phonon spectra based upon high angle measurements. 

Previous analysis of spin waves in LSFO [6] showed that the spectrum can be modeled ad- 
equately for energy transfers above 40 meV where phonon corrections are modest. We base 
our model on the facts that there are two different kinds of Fe ions, Fe^^ and Fe^+, and two 
nearest-neighbor exchange interactions which obey the Goodenough-Kanamori rules for the 
sign of the exchange; antiferromagnetism between half-filled Fe'^^-Fe^"'" pairs and ferromag- 
netism between half-filled and empty Cg orbitals in Fe'^'^-Fe^"'' pairs. The AW Heisenberg model 
Hamiltonian is 

H = -Jaf 2_^ Sj"*" ■ Sj"^ - Jf 2_^ Sj"^ ■ Sj"^ , (11) 

where sums are over each pair-type, S / and Sj represent the spin vector of the ?th and jth iron 
atom of the type indicated, and sums are over nearest-neighbor pairs with pairwise exchange 
values {J p or J af) determined by the charge ordered structure. 

Because of the small charge-transfer gap in 7?SF0, significant hybridization exists between 
Fe and O resulting in some fractions of doped holes residing on the oxygen site. Even con- 
sidering the presence of doped holes, the exchange between Fe'^+ and nominal Fe^"*" ions re- 
mains ferromagnetic. In the limit where the holes are on the iron site (corresponding to full 
Fe^"^ valence), F exchange occurs between half-filled and empty Cg orbitals according to the 
Goodenough-Kanamori rules. In the limit where a single hole is on the oxygen site, sharing of 
the spin-polarized oxygen electron also leads to F exchange. For the same reason, the presence 
of oxygen holes between Fe^"^ - Fe^^ pairs will reduce J af as compared to the parent insulator 
RFO. Furthermore, the spin values of Fe ions are affected by hybridization: if the hole is on the 
oxygen ion, the Fe oxidation state is lower and the effective spin of the Fe ions would be larger. 

Figure 7 shows the fitting results of PSFO and NSFO to the Heisenberg model of Eq. (11). 
Unlike the simple case of 7?F0, the fitting of AW Heisenberg model calculations and CEF ex- 
citations in i?SFO do not show quantitative agreement with the data especially at low energies 
where phonon corrections are more important. However, as explained in [6] the critical ratio 
\JfIJaf\ can be estimated by the splitting between the upper (> 70 meV) and lower energy (< 60 
meV) bands. For LSFO and PSFO, the WA^ Heisenberg model calculations show that these two 
main spin wave bands originate from F-like/AF-like spin waves that propagate along/between 
the metal centered domain wall, respectively. A rough estimate of the energy scale for these 
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FIG. 7: (color online) Comparison of /?SFO (R = Pr (a) and Nd (b)) magnetic scattering spectra at T = 10 
K (dots) and summed for 29 3°-30°. The Heisenberg model calculation (blue line), the CEF excitations 
(green line) and the sum of the two calculations (red line) are shown. 

excitations is, 



:5+\ 



E 



AF 



Ep = 3JFi2S'^ + S^ 
= 3 Jaf iJ + oJpS 






(12) 



Based on the previous reports on spin magnetic moments for different Fe ions, S^"*" is very 
stable and keeps as the value 5/2, while S^+ is likely to be larger than the atomic limit of 3/2 
due to hybridization and can be adjusted. [22, 23] The final fitting magnetic exchange energies 
are shown in Table VIII. 

The measured magnetic exchange ratio of Jp and | Jaf I can be compared to the theoretical 
predictions of the magnetic exchange mechanism for the CO in RSFO (R = La, Pr, and Nd) as 
shown in Fig. 8. In the magnetic exchange only model proposed by T. Mizokawa, et al.,[7] two 
charge ordered patterns are considered with alternating charge and spin ordering along either 
the (111) or (100) directions. Those with the observed pattern of charge ordering along the 
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TABLE VIII: Parameters of the AW Heisenberg model for RSFO {R = La, Pr and Nd) and the size of the 
charge transfer gap. [5, 24] Parameters of LSFO are from Ref. [6]. 





Lai/aSi-a/gFeOs 


Pi"i/3Sr2/3Fe03 


Ndi/gSra/gFeOs 


Spin momenta 








53+ 


2.5 


2.5 


2.5 


55+ 


2.0 ±0.1 


2.0 ±0.1 


2.0 ±0.1 


Magnetic exchange energies 








Jaf (meV) 


-3.5 ± 0.2 


-3.5 ± 0.2 


-5.5 ± 0.3 


Jf (meV) 


5.1 ±0.1 


5.1 ±0.1 


5.1 ±0.1 


ratio {\Jf/Jaf\) 


1.46 ± 0.05 


1.46 ±0.05 


0.93 ± 0.03 


chai-ge transfer gap 








A (meV) 


62 


58 


85 



cubic (lll)-direction are stable for large values of the exchange ratio, \Jf/Jaf\ > 1^ although 
the boundary between the two phases at T = 10 K depends on the value of NN exchange between 
Pe5+ _ pg5+ JQj^g i^y^^-) j-j^^j- jg present only in the (hypothetical) (100) charge ordered structure. 
According to the Goodenough-Kanamori rules, it is expected that this exchange is weakly AF 
due to the vr - bonding of half-filled t2g orbitals. In the limit where 755 ~ 0, the (1 11) order is 
stable when |y^ / Jaf\ > 1/2. A much more conservative estimate of 755 ~ | Jaf I results in the 
condition |7^ / Jaf\ > 1 for the stability of (1 1 1) charge order. LSFO and PSFO have exchange 
ratios that clearly favor the (111) ordering, even in the most conservative estimate for the value 
of 755. On the other hand, the exchange ratio for NSFO is slightly less than one. This suggests 
that the (1 1 l)-type charge ordering is less stable in NSFO as compared to LSFO or PSFO. This 
is consistent with the suppression of Tjy in the RSFO compounds. 

The reduction of the exchange ratio in NSFO largely arises from a ~ 20% increase of \Jaf\ 
as compared to LSFO and PSFO. The AF exchange in NSFO is comparable to that of the parent 
RFO compounds. Based on Table III and VIII, the increase of \Jaf\ in NSFO as compared to 
LSFO and PSFO could arise from differences in the CT gap and/or the lattice distortion of these 
compounds. The CT gap of NSFO is larger than those of LSFO or PSFO, hence electrons will be 
more localized and the magnitude of AF exchange energy of NSFO is expected to increase. In 
addition, the effect of lattice distortion on the magnetic exchange should also be considered. The 
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FIG. 8: (color online) The charged order phase diagram for /?SFO (/?=La, Pr, and Nd) as a function of 
nearest-neighbor magnetic exchange. The (111) and (100) types of charge order are illustrated respec- 
tivley in the coiTcsponding portions of the phase diagram. 

effect of the lattice distortions in the RVO parent compounds is well understood. The tolerance 
factor, whose deviation from one indicates the propensity for lattice distortion, decreases from 
LFO to NFO (see Table II). The larger lattice distortion in NFO results in smaller Fe-O-Fe 
bond-angles that weakens the AF superexchange. The tolerance factor of the i?SFO compounds 
is closer to 1 than the parent RFO compound and varies only weakly throughout the 7?SF0 
series. Therefore, the large change in J af is unlikely to arise from an average change of the 
doping induced structural distortions. 

V) Conclusion 



Using inelastic neutron scattering, we determined that the similar spin wave spectra of LSFO 
and PSFO consist of two energy bands separated by a large energy gap, while the two bands 
merge into one in NSFO. The full magnetic bandwidth is determined mainly by the ferromag- 
netic exchange energy, J p, between Fe^+ and Fe^''" ions, and is found to be similar for the 
different 7?SF0 compounds. The AF exchange energies between Fe^+ ions, \J af\, which con- 
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trols the splitting of the upper and lower magnetic bands are more sensitive to R substitution. We 
determine Jp and \Jaf\ by comparison to a Heisenberg model. The ratio of these exchanges is 
an indicator of the role that magnetism plays in the formation of the charge ordered state. While 
LSFO and PSFO are in the regime where magnetic exchange can stabilize the charge ordered 
state, the case for NSFO is not as clear. The much lower exchange ratio in NSFO may come 
from the increase of the charge transfer gap that is caused by the smaller Nd ion. 
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